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^ ■ Abstract 



In order to determine the Wigner function uniquely, we introduce a new 
condition which ensures that the Wigner function has correct marginal distri- 



> 

o 

l/^ ■ butions along tilted lines. For a system in N dimensional Hilbert space, whose 

o . 

^ ■ "phase space" is a lattice with N"^ sites, we get different results depending on 

whether N is odd or even. Under the new condition, the Wigner function is 
• determined if N is an odd number, but it does not exist if is even. 

^3 : I. INTRODUCTION 

I In the quantum theory, we can define the "weight" , although this weight may be negative, 

I and the expectation value of an observable is to be expressed as the average of a real function 

on phase space corresponding to the observable. This weight is well-known as the Wigner 

function |jl| and usually assumed to satisfy the conditions: 

(a) We can obtain the marginal distributions along the coordinate and momentum axes by 
integration. 

(b) The Wigner function is real valued. 

(c) The Wigner function includes the same information as the density matrix. 
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For a continuous quantum system with the coordinate Q and the momentum P, these 
conditions are described in the following forms: 

(a)The integration of the Wigner function W{q,p) over the variable p gives the expectation 
value of the density matrix p for the eigenvector \q) of Q 

W{q,p)dp = {q\p\q). 

Similarly, the expectation value of the density matrix p for the eigenvector \p) of P is 
obtained by integration of the Wigner function over the variable q 



W{q,p)dq = {p\p\p). 

(b) The Wigner function is real valued, 

W*{q,p) = W{q,p). 

(c) Usually this condition is written in terms of the Fano operator which is defined by 

W{q,p) = TT[A{q,p)p]. (1.1) 

Moreover it is imposed that the Hermite conjugate operator A^{q,p) is employed for the 
inverse relation 



dqdpA^{q, p)W{q, p) , 

-oo 



:i.2) 



although these equations are not equivalent to condition (c) but sufficient conditions for it. 
It is not difficult to show the orthogonality and completeness of the Fano operator. 
We can check that these conditions are satisfied by the function W{q,p) 



W{q,p) 
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For this Wigner function the Fano operator becomes 

1 



A{q,p) 



27rh 



dr 



ipr/h 



q + 



:i.4) 



Inversely, when we attempt to determine the Wigner function under these conditions we 
find ambiguity as is pointed out by many physicists 0- and in a previous paper A new 
condition is needed for the determination of the Wigner function. Bertrand and Bertrand 
essentially required correct marginals along tilted lines in addition to basic marginals 
{i.e. condition (a)) and showed that there exists only one solution Eq.( |1.3D . 

This condition is important not only for the determination of the Wigner function but 
also for quantum tomography. When we can get the marginal distributions on all tilted lines 
from observations, using the inverse Radon transformation we find out the Wigner function 
which has the same information as the density matrix. Specially, if we consider a pure state, 
we can obtain the wave function up to a phase factor from the observations. 

Now we consider a quantum system in a vector space spanned by eigenvectors | g) {q = 
0,1 ■ ■ ■ , N — 1) for the "coordinate" operator. It is assumed that these eigenvectors satisfy 
periodical boundary conditions 

\q + N)=\q). (1.5) 
The eigenstates of the "momentum" operator are defined by 

b) = ^ E ^"''"l^), (1.6) 



and the phase space of this system is a lattice with A^^ sites. Here, a; is a primitive Nth 
root of unity 

u = e N . 

When the conditions corresponding to conditions (a) ~ (c) are imposed, Cohendet et al. 
P gave the Wigner function for the case where is odd. For the case with A^=even 
number, we presented the Fano operator in the previous paper 0. These solutions are not 
unique and we showed that the general solution under the three conditions is obtained by 
multiplying an (A — 1)^ x (A — 1)^ orthogonal matrix to a special solution 0]. Leonhardt 
permitted arguments of the Fano operator to take values of half-integer besides integer 



values and found the Fano operator for A^=even number which satisfies the three conditions 
corresponding to conditions (a) ~ (c) |T0[. However, we assume the arguments of the Fano 
operator to take only integer, namely, eigenvalues of coordinate and momentum operators 
in this note. Wootters discussed the Wigner function for the case where N is prime and 
found the function which satisfies the property similar to the marginalization on tilted lines 



for continuous systems ITT 



In this note, first, for a continuous quantum system we propose a condition determining 
the Wigner function uniquely which is essentially equivalent to that of the Bertrand and 
Bertrand ^ but which is in a simpler form than they imposed. We easily extend this 
condition to the system with a lattice phase space. In both cases, we can see that the 
Wigner function under the new condition has correct marginal distributions along tilted 
lines in phase space. 



II. THE WIGNER FUNCTION IN A CONTINUOUS QUANTUM SYSTEM 

A. Behavior of the Fano operator under the rotation 

In this section we discuss the Wigner function for a continuous quantum system paying 
attention to the marginal distribution along tilted lines. In the classical theory, we can 
regard any tilted line in phase space {q,p) as the axis for new variable q' defined by 

q' = q cos 6 — p sin 9 
p' = qsm9 + p cos 9 

Given the distribution function p{q,p) for old variables, the distribution p'{q',p') for new 
variables is obtained from the simple transformation rule 

p'{q',p') = p{q,p) = p{q' cos 9 + p' sm9, —q' sm9 + p' cos9). (2.1) 

In the quantum theory, the transformation corresponding to the above linear transformation 
is generated by the unitary operator Rg 
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R, = e-^",(Q'+P'y\ (2.2) 

because we have 

Q ^ Q' = RgQRg^ = Q cose - P sine, (2.3) 

P ^ P' = RePRs^ = P cose + Q sine. (2.4) 

Under this transformation, the Fano operator A{q,p) is transformed to A'{q,p) following 
the same rule as usual: 

Aiq,p) A'{q,p) = ReA{q,p)R^\ 

From Eg. ( P-.l[ ), the Wigner function W'{q',p') of the same density matrix p is given by 

W\q',p') = TT[pReA{q',p')R^']. 

Since we consider the Wigner function as the distribution function in the classical theory, 
we require that the Wigner function has the same property as Eq.( p.l|) : 

W'{q'^p) = W{q' cose + p' sine, —q' sine + p' cose). 

Thus we have the new condition 

RgA{q,p)RQ^ = A(q'cos6' +psin6', —(/sin 6* +pcos6'). (2.5) 

We can prove that the Fano operator ( |1.4| ) satisfies this condition by using the Fano 
operator in the form 

^ roc 

A(g, p) = J ^ d3d^e^{*(^-^)-('3-^)>/^ (2.6) 

Inversely, can we determine the Winger function namely the Fano operator from these 
four conditions (a) ~ (c) and Eq.( p.5|) ? In the following subsections, we discuss this point. 
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B. the Fano operator under the new condition 



We start with rewriting conditions (a) ~ (c) in terms of the Fano operator. Using the 
relations Eq. (|l.l| ) and Eg .( 11.21) between the Wigner function and the Fano operator we get 
the following equations: 



-oo 
oo 



A{q,p)dp=\q){ql (2.7) 
k{q,p)dq=\p){pl (2.8) 
(b): At(g,p) = A(g,p), (2.9) 

/oo „ ^ 

(51 I At(g,p) I q^){q[ I A(g,p) I q'^)dqdp = 5{qi - q'^q'i ' ^2), (2.10) 
-00 

Tr[Aiq,p)A\q\p')]=Siq-q')Sip-p'). (2.11) 

Because of completeness of the operators e*®^/^ and e~^'^^^^ , we can expand the Fano 

operator with these operators, 

1 -I ~ I 

A{q,p) = —- / dQdVa{q,p; Q, P)e*SP/7ig-*Q7'A^ ^2.12) 

This expansion makes it easy to calculate the transformed Fano operator by the unitary 
operator Rg. In terms of the Fourier transformed coefficient a{s,t; Q,V) which are defined 

by 

~ais,t; Q^V) = ^ [e-^^^/^^P*/^a(g,p; Q,V)dqdp] , (2.13) 
three conditions (a) ~ (c) are described in the simple forms 

(a) : fi(s,0;Q,P) = 5(Q)5(P-s), (2.14) 

a(0,t;Q,P) =5(t- Q)(5(P), (2.15) 

(b) : ~a{s,t-Q,Vy = e-'^'''%{-s,-t--Q,-V), (2.16) 

/oo 
dsdt [5(s,t; Q,Vy~a{s,t- Q' ,V')] 
-00 

= S{Q-Q')6{V-V'), 

/oo 
dQdV [d*{s,t- Q,V)~a{s\t'-Q,V)] 
-00 

= 5{s- s')5{t-t'). (2.17) 



Since the operators e ^Q^/'* and e^^^f^ are changed to 

^^g«QP/7ij^-l _ gicos6»sin6»Q2/2giQsm6IQ/;igiQcoseP/ft ^2 Xg) 

under the transformation by Rq, we get the condition for a{s, t; Q, V) from the new condition 
Eq.(iD, 

= a(s cos 6* + t sin 6^, f cos 6* — s sin ; Q cos 6' — P sin ^, Q sin 6' + P cos 6) 

^ ^-iiQ cos e-V sin 9){Q sin 9+V cos 9) /2n j^2 20) 

For each point {s,t), we can choose the angle 6 such that scos^ + tsin^ = 0. 

cos6' 



sin 9 - 



X exp 



Then, from Eg. (|2.2Cl| ), we obtain 

tQ + sV -sQ + tV \ 
{tQ + sV)i-sQ + tPy 

Using condition (a) and the property of the 5-function, we can rewrite the right hand side 
of the above equation in a simple form and obtain 

which is equivalent to the coefficient for the original Fano operator given by ( |1.4[ ). 

Thus, under the new condition Eq.( p.5| ) we can determine the Fano operator uniquely. 

In this section, we have investigated the behavior of the Fano operator under the rotations 
induced by the unitary operator Rq in phase space. We have to extend the rotations to 
more general transformations in order to apply this method to the Wigner function on a 



lattice, because it is impossible to find a rotation of general angles which have one-to-one 
correspondence on lattice sites. Instead of the unitary operator when we apply the 
unitary operator Ra/s 

we have the transformation from the operators Q and P to the operator Qg and Pg, 

Q^Qg = RapQKl = AQ + BP, 

P^Pg = RapPKl = dP + cq, 

so that we get 

i^^^je-^^^/'^i?-^ = e^ABPV2^-^AQr/n^-^BPP/n^ ^2.21) 

R^f^e'^^/^R-;, = e'CDQ^I\''^''^l'^e'^''^l'^, (2.22) 



where coefficients A^B^C and D are 



A = cos(2ya^), 

B = -W-sin(2y^), 
V a * 

fey I — 
C = J-smi2^a/3), 

D = cos(2ya^). 



After the almost same calculations as the above, we can show that the Fano operator in 
Eq.( |1.4|) or Eq. (|2.6| ) satisfies the similar relation 

R^pA{q,p)R-'^ = A{Dq - Bp, -Cq + Dp). (2.23) 

If we can find the linear transformation with one-to-one correspondence among lattice sites, 
we can do the same thing on a lattice. 

In the next section we try to apply this method to finding the Fano operator on a lattice. 
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III. THE WIGNER FUNCTION ON A LATTICE 



For a system with a lattice phase space explained in Introduction, we play the same 
game as we did in the preceding section. 

We impose the following conditions (a)' ~ (c)' on the Fano operator A{q,p) {q,p G Z^) 

(a) ': Y: A(g,p) = k)(g|, (3.1) 

J2 Hq,p) = \p){pI (3.2) 

(b) ': A{q,py = A{q,p), (3.3) 

(c) ' : Y: {qi\Hq,Py\<l2){q[\A{q,p)\q',) = ^C^C' 



n r))/\(n' n')''] = 

N 



Tr[A(g,p)A(g',p')1 = 4C^4y' (3-4) 



which correspond to conditions (a) ~ (c) on the Fano operator in a continuous quantum 
system. Here, 6^.^} is the Kronecker delta on Z^, 



1 {q = q' mod A^) 
{q^q' mod A^) 



Owing to the expansion of the Fano operator by the operators e^^^l^ and e"*'^^''^, we 
wrote some properties in simple forms in the preceding section for a continuous quantum sys- 
tem. We introduce the matrices P and S corresponding to these operators. The eigenvectors 
for the matrix P are | (/) (g = 0, 1, ■ ■ ■ — 1), 

P\q)=uj'^\q), (3.5) 

and the eigenvectors for the matrix 5" is | p) (p = 0, 1, ■ ■ ■ — 1) which is defined by Eq.( |1.6|) 

S\p) =uj-'''\p). (3.6) 

In the representation where the matrix P is diagonalized, we can write the matrices P and 
5* explicitly 
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1 



: io : 

: ■■• 



1 ••• 0^ 



(3.7) 



(3.8) 

1 

1 

Using this representation, we can check that these matrices satisfy the commutation property 



V 



SP = ioPS, 



(3.9) 



which is similar to the commutation property 



We can consider matrices and on the lattice as the operators eT'^Q'^l^ and e'^^/'^, 
respectively. We expect that the same arguments as in the preceding section holds in this 
case. We expand the Fano operator in power series of matrices S and P; 



A(g,p)= E a(g,p;n,m)5'^P"^. 
Using these coefficients a{q,p;n,m), conditions (a)' ~ (c)' become 
(a)': E «(9'^>;^>"^) = ^'^"""C' 



1 



q&Zff 

(b) ' : a{q, p; n, m) — uj~'^'^a* {q, p; N — n, N — m) , 

(c) ': E (^*il'P'^n^^)^i(l'P'^k,l) = ^dl^^^Sl^J, 



q,peZ]si 



(3.10) 



(3.11) 

(3.12) 
(3.13) 

(3.14) 
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Introducing the Fourier transformed coefficient a{s,t;n,m) 



d{s,t;n,m) = ^ ^ u^'u P^a{q,p;n,m), (3.15) 



we can rewrite the above conditions in the simple forms: 

(a) ': a{s,0;n,m) = l^sZ^6lJ!}, (3.16) 

a(0, t; n, m) = ^^i^d^i^^ (3.17) 

(b) ': a{s,t;n,m) =LU~''"'d.*{N - s,N -t;N -n,N -m), (3.18) 

(c) ' : «*(^' m)d{s, t; k, I) = -^Si^Jsl^^J, 

y: ~a*{s,t-k,ira{s',t'-k,i) = ^jiydi^,\ (3.19) 

In the previous paper we showed that there exist many matrices satisfying conditions 
(a)' ~ (c)' and those matrices are transformed to each other by an orthogonal matrix. 

A. The new condition for the Fano operator 

Now we try to determine the Fano operator uniquely by imposing a similar condition to 
the condition Eq.(|Z^). 

Instead of the unitary operator Rg in Eq.( [2.23| ) we use the unitary matrix U which 
induces the transformation for matrices P and S 

P^UPU-^ = apS^P", (3.20) 
S^USU-^ = UsS^'P^, (3.21) 

where integers /i, z/, k, A and phase factors determined such that the matrices 

UPU^^ and USU~^ satisfy the same conditions as matrices S and P do. From Eq.( p.9|) , for 
integers /t and A, we get 

Kz/-^A = 1, (3.22) 
and from 5*^ = 1 and P^ = 1, we have 
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,(Af-l)W2 



{N~1)kX/2 



(3.23) 

(3.24) 



Comparing Eqs.( 3.20 ) and ( p.21 ) with Eqs.( 2.21 ) and ( 2.22 ), we find the correspondence 
for the parameters of both transformations 

A ^ u, B ^ —/X, 
C ^ -\ D^K. 

Thus, on the Fano operator on the lattice, we impose the new condition 



UA{q,p)U ^ = A{K,q + iip,up + \q). 



(3.25) 



B. the Fano operator under the new condition 



First, we find the condition for the coefficient a{q,p] n, m) in Eq. (|3.10| ) which is equivalent 



to Eq.(^. 

From Eq. (|3l0|) , the LHS of Eq.(^) becomes 



UA{q,p)U-^= J2 a{q,p;n,m){USU-^Y{UPU-^T 
= a{q,p]n,m)ay;;{S''P^)''{S^'P'')'^ 
= a(g,p;ri,m)u;'^("'™)5''"+^™P^"+^'", 



n,m£Zjsi 



where (p{n, m) is given by 



(3.26) 



(f){n, m) = -{n{N — n)«;A + m{N — m)fii/} — nmfiX, 
and we can show that periodic function 

^Hn+aN,m+bN) ^ ^</.{n,m) (^^^ ^ ^ integer). 

Instead of the summation with respect to n and m, we make the summation with respect 
to n' and m' which are given by 
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n' = nn + fim 
m' = Xn + vm 



(3.27) 



Owing to the condition Eg. ( p.9|) , we can solve Eq.( p.27|) for n and m, 



n = un' — fim' 



(3.28) 



m = —Xn' + Km' 
Substituting these equations into 0(n, m) we obtain 

(f){iyn' - fim', -Xn' + Km') = (p'{n', m') + 6 + N x (integer). (3.29) 

where (j)'{n',m') and 6 are given by 

(f)'[n', m') = ^{uXn'{N - n') + fiKm'{N - m')} + fiXn'm', (3.30) 
5 = AT I ^^^iHzJi^n' + "^^"-^-^^mj . (3.31) 

We can see that the values i'X{k — /i — 1) and Kfi{u — A — 1) are even under the condition 
Eq.( |3.22 ) as follows. We have two cases where i>X{k — /i — 1) is odd. The one is the case 



where k,X, fi and u are odd, the other is the case where A and u are odd and k and ^ are 
even. For both cases the value ku — fiX becomes even and the condition Eq. ( p.22| ) is not 
satisfied. Similarly Kfi{u — A — 1) in the second term is even, so that we can replace the 
function cj<^('^"'-M'»'--An'+«m') by u;<^'("''™'). It can be easily checked that the function cj<^' 
is a periodic function of period A^, 

^^'in'+aN,m'+bN) ^ (^^ ^ = integer), 

and that this function is a function on Zn x Zn- As the coefficients of the linear transfor- 
mation Eq. ( |3.27D and the inverse transformation Eq. ( |3.28| ) are integer, the map defined by 



Eq.( |3.27D has one-to-one correspondence on Zn x Z^. So we can make summation by n' 
and m' over Z^r x Z]^ instead of n and m over the same region: 

UA{q,p)U-' = a{q,p;iyn' - iim',-Xn' + Km')cu^'^'''^'^'^S'''P"'' (3.32) 
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Thus, the condition (|3.25|) in terms of a{q,p; n, m) is written by 

a{q,p; un — fim, —Xn + Km)uj'^ '•"''^^ = a{Kq + fip, vp + Ag; n, m). (3.33) 
For the Fourier transformed coefficient a(s, t; n, m), this condition becomes 

a{ys + At, /iS + nt; n, m) = u'^ ^"''''^^a{s, t; un — fim, —Xn + nm). (3.34) 
Replacing variables s and t hj ks — Xt and ut — fis, respectively, we get 

a{s, t; n, m) = u'^ ^"^'"^^a^KS — Xt, ut — /is, t; un — fim, —Xn + Km). (3.35) 

As the coefficient a{s,t;n,m) is a function on Zj\f x Z^r x Z^v x Zj\f, if we can choose the 
integers k and A such that 

Ks — Xt = integer x A^, (3.36) 

owing to condition (a)', we can obtain a{s,t]n,m). Thus, we can determine the coefficient 
a(s, t; n, m) for any integers s and t, uniquely. 

Now, we try to get an explicit form of a{s,t;n,m). We choose a and r as A and k 
respectively, which are defined by. 



s = ^a + N 



s 

N 
t ■ 

N. 



(0<a<A-l), (3.37) 
(0<r<iV-l), (3.38) 



where ^ is the greatest common divisor of and f r, namely, a and r are relatively prime 
to each other. If we choose a and r as k and A, respectively, from number theory, it is clear 
that there exist the integers /i and u satisfying the equation ( |3.22| ). Then we obtain 

a{s, t; n, m) = u'^ '•"■''"•*a(0, ^; vn — /im, —Xn + Km) 

= j-y^-'Si^C- (3-39) 
From Eg. ( p.30| ), the (f)'{s,t) with the choice of k and A becomes 
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t) = ^{maT{iy + ^ + 1) - NCar - ^Vr}. 

Using the same discussion as we did before Eq. p.32| ), we can show that the first term in the 
above equation is multiplied by an integer, and we have 

^4''is,t) ^^-limar+e^r)^ (3.40) 

We can replace ^ by in the first term in the phase factor u'^'^^'''^\ since uj^^'^ = —1 and 
0^^5/2 jg gq^g^i a;^«'/2. We get 

~a{s, t; n, m) = (3.41) 



On the process to get this solution, we imposed the condition ( ^.341 ) for the integers k. A, 



/i and 1/ corresponding to the only one line which passes through the fixed point (s, t). So it is 
not clear that this solution satisfies the condition (|3.34|) for the four integers corresponding 



to other lines. We have to investigate whether the four conditions are satisfied by this 
solution or not. 

When is odd, we can drop the term proportional to in the phase factor as A^ + 1 is 
even, and the solution becomes simple 

which is equivalent to the solution given by Cohendet et al. ||^ 

a(,M;«,m) = l = (3.43) 

I ^.-^Sim^ (n = odd) 

It is proved that this solution satisfies conditions (a)' ~ (c)' and we can easily check that 

the coefficient d{s,t;n,m) satisfies the condition Eq. (|3.34|) . 

However, when A^ is even, the solution Eq. (|3.4l|) does not satisfy the condition (b)' 

because 



7^- 



S 

N 



N 

Thus there dose not exist the Fano operator which satisfies original three conditions (a)' 
(c)' and the new condition Eq. (|3.25|) . 
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IV. SUMMARY AND DISCUSSION 



In this note, in order to determine the Wigner function uniquely, we proposed a new 
condition. For a continuous quantum system, the condition Eq.( p.23| ) is the almost same 
condition as Bertrand and Bertrand presented [Q. However, this condition is simpler than 
their condition and can be easily extended to the system with a lattice phase space. 

If the Fano operator satisfies the condition Eq.( p^.23 ), we can easily estimate the integra- 



tion along a tilted line, say Cq + Dp = po, and we have 

A(g = Dr — Bpo,p = Apo — Cr)dr 

-oo 

Raf3A{po,r)R~^dr 

-oo 

/" oo „ 
A{po,r)dr 
-oo 

= Rap I Po){PO I Kp- 



R-al 



Clearly, the state i?Q,/3 | p) is an eigenstate of the Hermite operator Fg, 

Pg = RafsPR-p = CQ + DP. 

The integration of the Wigner function along this tilted line gives the distribution function 
over the eigenvalues of the operator CQ + DP. Thus the Wigner function gives a cor- 
rect marginal distribution along a tilted line and this is important property for quantum 
tomography. 

On the analogy of a continuous quantum system, we presented the condition Eq.( p.25D 



for the system whose phase space is a lattice with A^^ sites. Under this condition, for the 
case where is odd, we found out the same solution that Cohendet et al. did 0, and for 
the case where is even we showed that there does not exist the solution. 

For the Fano operator on a lattice which satisfies the condition Eq. (|3.25|) , the summation 
of the Fano operator over the sites on the line np — Xq = po becomes 

J2 Ml = K,r + fipo,p = upo + Xr) = U \ Po){po I U~^, 
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which is the projection operator to an eigenvector of USU~^ as is seen from Eqs.( p.6|) and 
( p.21| ), so that we can get the correct marginal distribution on a tilted line for this system. 

Finally, we point out that the Fano operator is determined uniquely from only two con- 
ditions, namely, condition (a) and Eq. ( p.23| ) for a continuous quantum system and condition 



(a)' and Eq.(|3.25|) for a lattice. 
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